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Abstract 

We have recalculated the interior structure of very massive stars of uniform chemical composition 
with the OPAL opacity. Very massive stars are found to develop a core-halo structure with an extended 
radiative-envelope. With the core-halo structure, because a more massive star has a more extended 
envelope, the track of the upper zero-age main-sequence (ZAMS) curves redward in the H— R diagram at 
> 100Af o (Z = 0.02), > 7OAf (Z = 0.05), and > 15M Q for helium ZAMS (X = 0., Z=0.02). Therefore, 
the effective temperatures of very massive ZAMS stars are rather low: e.g., for a 2OOAf star, logT c ff = 4.75 
[Z = 0.004), 4.60 (Z = 0.02), 4.46 (Z = 0.05), and 4.32 (Z = 0.10). The effective temperatures of very 
luminous stars (> 12OAf ) found in the LMC, the SMC, and the Galaxy are discussed in relation to this 
metal dependence of a curving upper main-sequence. 

Key words: Stars: abundances — Stars: individual(the Pistol star) — Stars: interiors — Stars: 
massive — Stars: supergiants 



Introduction 

Recent infrared observations have found a number of 
very luminous stars in young clusters near to the galac- 
tic center, in 30 Doradus in the LMC and in star-forming 
regions in the SMC (Nagata et al. 1993, 1995; Najarro 
et al. 1997; Messey, Hunter 1998). The Pistol star, a 
member of the AFGL 2004 young cluster (Nagata et al. 
1993; Figer et al. 1996), is one of the brightest stars 
known in the local group of galaxies. The luminosity 
and the temperature of the Pistol star were derived to 
be logL(L ) = 6.6 and log T e g(K) = 4.15, respectively, 
from a near-infrared data analysis (Figer et al. 1998). 
Very luminous stars with such a low surface tempera- 
ture are rarely found in conventional observations, which 
stimulates a theoretical interest in its evolutionary stage. 
Compared with the evolutionary path of very massive 
stars in the H-R diagram, Figer et al. have derived the 
initial mass of the Pistol star to be 2OO-25OAf and the 
age to be 1.7-2.1 Myr. In the LMC, there also found 
a number of very luminous 03 stars in the R136 clus- 
ter in 30 Dor, several of which have a mass in excess of 
120 M Q (Messey, Hunter 1998). One of them, Melnick 
42, is analyzed in detail by Pauldrach et al. (1994), who 
estimate the luminosity and the mass of the star to be 
logL(L ) = 6.6 and M = 150M Q . 



In this way, very luminous stars are found in various 
chemical circumstances. Their inferred mass, age, and 
evolutionary status provide fundamental information for 
studies of star-forming regions. In order to estimate their 
age and mass, we require theoretical tracks of very mas- 
sive ( > 12OAf ) MS stars of both Populations I and II. 
With the discovery of the very massive stars (L > 1O 6 L ) 
mentioned above, we need to reanalyze the structure of 
very massive MS stars. 

Kato first found from an analysis that very massive 
stars develop a core-halo structure, which results in a 
redward bending of the upper main-sequence. Kato 
(1985) calculated the structures of very massive Newto- 
nian stars with Compton-scattering opacity, and found 
that vary massive stars develop an extremely extended 
radiative envelope. Such a structure is caused by an out- 
ward increase of the opacity in the radiative region. With 
the Compton-scattering opacity, this core-halo structure 
appears only in very massive stars (> 10 6 M Q ). When 
a star develops a core-halo structure, the photospheric 
temperature decreases owing to the extended envelope. 
Because a more massive star has a more extended en- 
velope, the upper part of the main-sequence curves to 
the right in the H-R diagram at very high luminosity, 
L > 10 10 L Q , for Pop I (Kato 1986). Such a massive 
star of ~ 1O 6 M may be a theoretical problem rather 
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than a realistic one; but after the new opacity, we ex- 
pect that the core-halo structure and the resultant main 
sequence bending will be realized in much less-massive 
stars. This is because the Compton-scattering opacity 
increases outward only a few percent in the radiative re- 
gion, whereas the new opacity has a prominent peak at 
around T~2x 10 5 K that must cause the structure to 
change effectively. 

The bending of the main-sequence appears for a much 
smaller mass in the case of helium stars with the Los 
Alamos radiative opacity. Langer (1989) has calcu- 
lated chemically uniform helium stars with central he- 
lium burning as models of Wolf-Rayet stars. The pure 
helium main-sequence (Y = 1.0) only indicates bend- 
ing at > 60M Q , but sequences of C/O-rich helium stars 
curve rightward in less-massive stars (e.g., at > 15M Q 
for Y = 0.02, C + O = 0.98). The bending appears in 
less-massive C/O-rich stars because of a larger radiative 
opacity of C/O-rich matter. 

The core-halo structure is also reported in helium 
main-sequence stars with optically thick winds. Kato 
and Iben (1992) have presented the interior structure 
of helium main sequence stars with artificially enhanced 
opacity in a model of Wolf-Rayet stars. The core-halo 
structure is developed in stars of mass 15 - 30 M Q , in 
which the radiative envelope is extended in a way that the 
density profile changes while responding to any opacity 
variation. In the H-R diagram, the main-sequence runs 
from the lower-left to the upper-right, contrary to the 
usual main-sequence. This main-sequence corresponds 
to the upper half part of the curved main-sequence de- 
scribed above. 

In this way, chemically homogeneous stars show a 
core-halo structure which results in the curved main- 
sequence when the opacity monotonically increases out- 
ward. The new opacity, which varies around a prominent 
peak, requires a recalculation of the core-halo structure 
and bending of the main-sequences. Therefore, we have 
recalculated the interior structure of very massive main 
sequence stars with the OPAL opacity. After the new 
opacity appeared, there have been presented many cal- 
culations on massive-star evolution for various problems, 
such as evolution with wind mass-loss, instability against 
radial oscillation, WR star modeling, evolutionary con- 
nection between O stars, LBV, and WR stars etc. (e.g., 
Schaller et al. 1992; Heger, Langer 1996; Glatzel, Kiri- 
akidis 1998; Stothers, Chin 1996, 1997, and references 
therein); but little attention has been paid to the core- 
halo structure and the resultant bending of the main- 
sequence. The main-sequence up to 120 Mq shows no in- 
dication to turn to the right for Z < 0.02, (Schaller et al. 
1992), and it slightly does so for Z = 0.03 (Glatzel, Kiri- 
akidis 1993); however, more massive stars up to 300M Q , 
ZAMS curves redward at > 120M Q for Z = 0.02 and 
0.04 (Figer et al. 1998). In the present paper, we have 



concentrated on massive zero-age main-sequence stars 
in order to examine the basic properties of the core-halo 
structure and to confirm the bending of MS. The next 
section describes the method and assumptions of the nu- 
merical calculation, and section 3 shows the core-halo 
structure of chemically homogeneous stars and presents 
curved ZAMS for Populations I and II in the H-R di- 
agram. A comparison with observational data of very 
luminous stars, such as the Pistol star, is given in discus- 
sions. 

2. Calculations 

The structure of massive stars with uniform chem- 
ical composition is obtained by solving the equations 
of hydrostatic balance, mass continuity, energy trans- 
fer by diffusion and by convection, and energy conser- 
vation with the assumption of spherical symmetry. We 
calculated two sets of ZAMS solutions, i.e., with/without 
mass-loss. The wind mass-loss is assumed to be in the 
quasi-steady state in which the heat flux is steady in 
the g-coordinate, i.e., the gravitational contraction term 
in the energy-conservation equation is approximated by 
c g = -T(M/M)(5s/8 In q) u where q = M r /M is the 
mass within radius r divided by the stellar mass, and 
the suffix t denotes the partial derivative with constant 
time. These equations and assumptions are essentially 
the same as those in Kato (1980). This quasi-steady 
state condition stands well for a chemically uniform star 
with mild mass-loss. We have checked it by reproducing 
the diffusive luminosity distribution obtained by a hy- 
drodynamic calculation for a 15 Mq He MS star (Heger, 
Langer 1996). The wind mass-loss rate is assumed to 
be zero and 5 x 10~ 5 Af Q yr~ 1 for MS stars. Because 
of the effects of wind mass-loss on the stellar structure 
is very small, the interior structure hardly changes up 
to 1 x 10 _4 M Q yr _1 : the effective temperature increases 
only by AlogT ff = 0.003 if we include a mass loss of 
M = 1 x 10- 4 M Q yr - 1 in the 200 M Q model. The up- 
dated OPAL opacity (Iglesias, Rogers 1996) is used and 
the mixing length parameter for convective energy trans- 
port is set to be 1.5. The chemical composition of stars 
is assumed to be uniform with X = 0.7 for hydrogen, 
and Z = 0.004, 0.02, 0.05, and 0.1 for heavy elements by 
weight. In addition to these compositions, we have calcu- 
lated additional models with different combination for a 
comparison: X = 0.35 and Z — 0.05, and helium ZAMS 
of X = 0, and Z = 0.004,0.02,0.05, and 0.1 without 
mass loss. 

3. Interior Structure of Massive Stars and the 
H— R Diagram 

Figure 1 shows the distributions of the density and the 
temperature of chemically homogeneous stars of 40 and 
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Fig. 1.. Temperature and density distributions of chemically 
uniform stars (X = 0.7, Z = 0.02) of 40 and 200 Mq (solid 
curve). The dot denotes the outer edge of the convective core. 
The outermost point of each curve denotes the photosphere. 
The dashed and dotted curves denote the stars of 200Mq with 
(X = 0.7, Z = 0.004) and (X = 0.7, Z = 0.10), respectively. 
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Fig. 2.. Variation in the diffusive luminosity L r , the local Edding- 
ton luminosity, Lsdd = 47rcGM r /K, and the opacity k for mod- 
els in figure 1, except the 2OOM with (X = 0.7, Z = 0.004), 
which is omitted from this figure. L r , solid curve; Z/Eddi dotted 
curve; ft, lower solid curve. 



200M© for Z = 0.02. The 2OOM star has a quite differ- 
ent structure from the 40 Mq, because it develops an ex- 
tended isothermal radiative-envelope where the density 
is almost constant. Such a core-halo structure develops 
well in very massive stars (> 15OM ), but does not do so 
in less-massive stars, such as the 4OM , as shown in this 
figure. These two different types of homogeneous stars 
have already been pointed out by Kato (1985). Follow- 
ing her way, we call the core-halo structure a Type-II 
solution, and for the other one, the usual structure like 
in a 40M© star, a Type-I solution. 

Figure 2 shows the distributions of the diffusive lumi- 



Fig. 3.. Same as those in figure 1, but for He main-sequence stars 
of 12 and 4OM with (X = 0., Z = 0.02). The thin dashed and 
dotted curves denote the stars of 40Mq with Z = 0.004 and 
Z = 0.05, respectively. 



nosity, the local Eddington luminosity, 

^Edd = Z , 



(1) 



and the opacity. As shown in this figure, the local Ed- 
dington luminosity decreases outward in the radiative re- 
gion, corresponding to an increase of the OPAL opacity 
toward the peak at T ~ 2 x 10 5 K. The 200M Q star 
shows the super-Eddington luminosity in the outer ra- 
diative region, where the convective heat transport is in- 
efficient and a wide isothermal region develops to form a 
core-halo structure, as shown in figure 1. In the 40M© 
model, the diffusive luminosity does not reach the Ed- 
dington luminosity in the radiative region, and thus no 
core-halo structure appears. 

We also calculated the structures of helium ZAMS 
stars (X = 0.0) and helium-rich stars (X = 0.35, 
Z = 0.05). Type-II solutions also appear in these stars, 
and their basic properties are the same as those in ZAMS 
stars. Figure 3 shows the density and the temperature 
distributions for helium stars of 12 and 40M Q ; the latter 
shows the core-halo structure. 

Such a core-halo structure is more extended in metal- 
rich stars. Figure 1 demonstrates the difference of inte- 
rior structures of a 200M Q star with various metal con- 
tents, Z = 0.004, 0.02, and 0.10. Metal-rich stars develop 
a wide isothermal radiative envelope. As the envelope is 
most extended, the effective temperature of the star is 
the lowest of these three stars. 

Table 1 gives the physical quantities of solutions for 
very massive stars with various sets of chemical compo- 
sition, X and Z. The four columns next to Z give the 
photospheric values of the radius R p h, temperature T p h, 
luminosity L p h, and ratio of the diffusive luminosity to 
the Eddington luminosity at the photosphere £ p h/£Edd- 
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Table 1. Characteristic values of very massive star. 
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lno-T (K\ 
10gi c (i^) 


40 


0.7 


0.02 


9.3 


4.62 


5.35 


0.65 


0.38 


7.57 


60 


0.7 


0.02 


12 


4.65 


5.71 


0.74 


0.28 


7.59 


100 


0.7 


0.02 


18 


4.67 


6.11 


0.84 


0.19 


7.61 


200 


0.7 


0.02 


40 


4.60 


6.57 


0.92 


0.091 


7.63 


250 


0.7 


0.02 


55 


4.57 


6.71 


0.95 


0.077 


7.63 


300 


0.7 


0.02 


71 


4.54 


6.82 


0.96 


0.066 


7.63 


500 


0.7 


0.02 


160 


4.44 


7.10 


0.97 


0.042 


7.65 


1000 


0.7 


0.02 


3000 


4.08 


7.46 


0.95 


0.029 


7.66 


100 


0.7 


0.05 


26 


4.58 


6.09 


0.87 


0.19 


7.58 


200 


0.7 


0.004 


20 


4.75 


6.58 


0.87 


0.096 


7.66 


200 


0.7 


0.05 


76 


4.46 


6.56 


0.93 


0.091 


7.60 


200 


0.7 


0.1 


140 


4.32 


6.55 


0.89 


0.10 


7.58 


40 


0.35 


0.05 


15 


4.59 


5.68 


0.85 


0.24 


7.58 


100 


0.35 


0.05 


75 


4.41 


6.34 


0.86 


0.10 


7.62 


150 


0.35 


0.05 


200 


4.26 


6.59 


0.89 


0.073 


7.63 


8 


0. 


0.02 


0.82 


5.04 


4.94 


0.73 


0.41 


8.13 


12 


0. 


0.02 


1.1 


5.07 


5.29 


0.88 


0.32 


8.28 


15 


0. 


0.02 


1.3 


5.08 


5.47 


0.93 


0.28 


8.30 


30 


0. 


0.02 


8.3 


4.79 


5.96 


0.94 


0.18 


8.31 


40 


0. 


0.02 


20 


4.64 


6.14 


0.97 


0.15 


8.32 


60 


0. 


0.02 


71 


4.43 


6.38 


0.80 


0.12 


8.33 


100 


0. 


0.02 


670 


4.02 


6.67 


0.02 


0.096 


8.36 


130 


0. 


0.02 


850 


4.00 


6.81 


0.017 


0.087 


8.36 


40 


0. 


0.004 


2.2 


5.13 


6.14 


0.99 


0.15 


8.32 


40 


0. 


0.05 


130 


4.25 


6.14 


0.81 


0.15 


8.32 



* The ratio of the photospheric luminosity to the Eddington luminosity at the surface, 
t The ratio of the mass of the radiative envelope to the total stellar mass. 
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Table 2. Critical mass between Type I and Type II solutions 
of ZAMS star. 



z 


Mass (M©) 


logT eff (/0 


log L pb (L 


0.004 


200 


4.75 


6.58 


0.02 


90 


4.67 


6.03 


0.05 


70 


4.60 


5.81 


0.1 


60 


4.54 


5.66 



Table 3. Critical mass between Type I and Type II solutions 
of Helium ZAMS stars. 




E-r 

o 



Mass (M ) log Tea (K) logL ph (L Q ) 



0.004 


30 


5.13 


5.96 


0.02 


15 


5.07 


5.29 


0.05 


10 


5.03 


5.14 


0.1 


7 


4.98 


4.83 



[-7/7] Soi 



The next gives the ratio of the mass of the radiative en- 
velope to the total stellar mass, 1 — (M c /M), where M c 
denotes the mass of the convective core. The last column 
gives the central temperature T c of the star. 

In the main-sequence with (X = 0.7, Z = 0.02) in ta- 
ble 1, T p h increases with the stellar mass until it has the 
maximum value at ~ IOOMq, and decreases after that 
as the core-halo structure develops. Here, we define the 
critical mass which divides the Type-I and Type-II solu- 
tions as the stellar mass of the maximum effective tem- 
perature in each sequence. This critical mass is tabulated 
in tables 2 and 3, for ZAMS and He-ZAMS, respectively. 

Figure 4 shows four tracks of ZAMS with X = 0.7, and 
Z = 0.004,0.02,0.05, and 0.10, left to right, in the H-R 
diagram. The metal-poor main-sequence (Z = 0.004) 
curves slightly rightward at > 300M Q , while the metal- 
rich sequences curve strongly at above several tens of M Q . 
This is because the core-halo structure develops well in 
metal-rich stars, and the effective temperature is lower 
than that of the metal-poor star with the same mass. 
Therefore, the effective temperature depends strongly on 
the metal contents, while the luminosity depends weakly 
on the metallicity, as shown in figure 4. 

Figure 5 depicts U—V curve of the stars of 40, 200, and 
as the extreme massive case 1000 M Q , to demonstrate the 
difference of in the interior structures of the Type-I and 
Type-II solutions, where U and V are the homologous 
invariants, defined by 



Fig. 4.. Track of zero— age main-sequence with various metal 
contents in the H— R diagram. The solid curves denotes the 
main-sequence with X = 0.7, and Z = 0.004,0.02,0.05, and 
0.10, from left to right. The stellar mass is noted next to the 
dot in units of Mq. The dotted curve up to 240Mq denotes 
ZAMS (Z = 0.03) with old opacity (Maeder 1980). 



and 



V = GM r p/(rP). 



(3) 



U = 4irr 3 p/M r 



(2) 



Here, U represents the density divided by the mean den- 
sity, and V is the ratio of the gravitational energy to the 
thermal energy. A Type-I solution, such as of the 40M Q 
main-sequence star, shows almost a monotonic increase 
of V from the center to the photosphere, whereas the 
Type-II solutions of 200 and 1000 M Q make a deep dip or 
a loop around V ~ 2. The decrease of V toward the dip is 
caused by a quick decrease in the density outward, which 
corresponds to a steep rise in the opacity. After the loop, 
the opacity begins to decrease, which keeps the density 
at a relatively large value, and then V increases again. 
Note that such a loop structure is a characteristic prop- 
erty of red giants with hydrogen-shell burning, which has 
an extended envelope around a dense core (Hayashi et al. 
1962). 

4. Discussion 

Figure 6 shows ZAMS tracks for various chemical com- 
position. The thick and thin curves denote the ZAMS 
with X = 0.7 and Z = 0.0, respectively. One additional 
sequence of X = 0.35 and Z = 0.05 is also shown for 
a comparison. Table 1 shows that a hydrogen-deficient 
star (X = 0.35, Z = 0.05) has a higher luminosity and a 
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Fig. 5.. U-V curves for main-sequence stars of 40, 200, and 
1000 M with X = 0.7 and Z = 0.02, where U = 4vrr 3 p/M r 
and V = GM r p/(rP). The stellar mass is attached beside the 
curve. Two models of 200 Mq with Z = 0.004 (dotted curve) 
and Z = 0.10 (dashed curve) are also shown. The outer edge of 
the convective core is indicated by filled circles for 1000, 200, 
and 40 Mq from left to right. 

lower temperature, compared with that of a hydrogen- 
rich star (X = 0.7, Z = 0.05) of the same mass. This 
makes a good contrast to the weak dependence of the 
luminosity on Z; i.e., the stellar luminosity is almost 
independent of Z for a given mass. The difference be- 
tween helium-rich stars (X = 0.35) and solar composi- 
tion stars (X = 0.7) in luminosity and temperature is 
AlogL = 0.22 and AT cff = -0.26, for 150 M Q , and 
AlogL = 0.195, and AT off = -0.50 for 200M Q . In other 
words, observational information on helium enrichment 
is important when we derive the stellar masses from the 
observed luminosities. 

We now compare our theoretical ZAMS tracks in the 
H-R diagram with several very luminous stars which 
were recently discovered. Figure 6 also depicts the po- 
sition of several very bright stars [log L/(Lq) > 6] in 
our Galaxy, the LMC, and the SMC. The distribution of 
the LMC and the SMC stars (filled symbols) is consis- 
tent with our ZAMS curve of Z = 0.04, because Magel- 
lanic stars are known to be metal poor. A very massive 
LMC star, Melnick 42, denoted by the filled circle, has 
been examined in detail by Pauldrach et al. (1994), who 
derived the stellar parameter to be log(L/L Q ) = 6.6, 
T cff = 50500 K, Z = Zq/4, and M/M = 150, from 
model fitting of non-LTE UV spectrum with HST UV 
spectrum. The position of this star in figure 6 (filled 
circle) is close to the ZAMS of Z = 0.004 and quite con- 
sistent with our aspect. From our sequence of X = 0.7 
and Z = 0.004, its mass is estimated to be ~ 200M Q ; if 
the star is helium-rich, the mass is slightly smaller than 
this value, as mentioned above. 

The distribution of these Magellanic Cloud stars seems 



Fig. 6.. Theoretical ZAMS (X = 0.7) plotted in the H-R 
diagram (thick curves). The heavy element content of each 
sequence is assumed to be Z = 0.004, 0.02, 0.05, and 0.10, from 
left to right. The two horizontal dashed curves connect the 
points of 200 Mq (above) and of 100 Mq on each ZAMS track. 
The uppermost point of each curve corresponds to 800, 1050, 
400, and 280 Mq. One additional track of helium-rich ZAMS 
(X = 0.35, Z = 0.05) for stars from 20 Mq to 200 M is 
denoted by the dot-dashed curve. Very luminous stars in the 
Galaxy (open symbols) and in the Magellanic Clouds (filled 
symbols) are also plotted. The two open circles with the error 
bar denote the two possible positions of the Pistol Star (Figcr ct 
al. 1998), the open circle r\ Car (Humphreys, Davidson 1994), 
and the filled circle Melnick 42 (Pauldrach et al. 1994). Several 
O stars in the Galaxy, the SMC, and the LMC (Puis et al. 
1996) arc plotted by the open squares, the filled triangles, and 
the filled squares, respectively. The four open triangles denote 
bright He-rich stars in the galactic center region (Najarro et 
al. 1997). 



to be on the whole slightly leftside to that of the galactic 
stars (open square). This is consistent with the metal de- 
ficiency of the Magellanic Clouds stars, although the data 
number is not sufficient for making any definite state- 
ment, and some of which possibly have left the main 
sequence to a redward evolutionary excursion. 

This figure also shows four massive stars in the galactic 
center (open triangles). Najarro et al. (1997) have ana- 
lyzed He I lines of these stars by a non-LTE radiative- 
transfer method for a pure H He spherical atmosphere. 
They showed that all of them are strongly He enriched 
(Hc/H > 1). Because these IRS stars are in the galactic 
center, they are possibly enriched in metal as well. In 
the H-R diagram, they are located close to our theoret- 
ical curve of X — 0.35 and Z = 0.05 in the upper-half 
part that represents core-halo structure. Therefore, if 
their interior structure is not far from that of our uniform 
model, we can expect that their low effective temperature 
may be attributed to the core-halo structure. 

Another candidate for the core-halo structure is the 
Pistol star. The luminosity and temperature of the Pis- 
tol star has been derived by Figer et al. (1998) from a 
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near-infrared observation to be two possibilities: (case 
L): log L(Lq) = 6.6 ± 0.2 and logT eff (i<r) = 4.15 ± 0.01 
K and (case H): logi(L Q ) = 7.2 ± 0.2 and logT off (K) = 
4.33 ± 0.01 K. They have analyzed the star to be slightly 
helium rich, nne/^H = 0.13-0.12; This gives X = 0.63- 
0.66 for Z = 0.02-0.05. The position of case L is very 
close to the curve of ZAMS with X = 0.7 and Z = 0.1. 
This metallicity seems to be too high, but we may not 
exclude Z = 0.1 for the Pistol star, because it is in the 
star-forming region in the galactic center. If the Pistol 
star is extremely metal-rich, the star can be interpreted 
to be a near zero-age main-sequence with a core-halo 
structure. In the case of Z = 0.02-0.05, the star is inter- 
polated to be a young star that is evolved redward apart 
from ZAMS. In any case, the stellar mass is estimated 
to be ~ 250M©, which is in a good agreement with 200- 
250 Mq obtained by Figer et al. In case H, the star is 
close to the main-sequence with Z = 0.02 and the mass 
is estimated to be extremely high, as much as 700 Mq. 

In this way, the low effective temperature of a very 
luminous star can be interpolated in part by the core- 
halo structure. It is difficult, however, to confirm our 
prediction, because we need to distinguish a very young 
star having a core-halo structure from an evolved star. 
Observational estimates of the helium and heavy element 
contents of such stars are useful information as well as a 
detailed evolutional calculation. 

Figer et al. (1998) have shown ZAMS tracks up to 
3OOM with Z = 0.02 and 0.04, of which the upper half 
curves rightward, which is qualitatively in good agree- 
ment with ours. Their temperatures are systematically 
higher than those of ours (AlogT cff ~ 0.1 at 2OOM ), 
whereas the luminosity is the same. This difference may 
be caused by a different treatment of the radiative en- 
velope, because they combined an interior solution to a 
steady wind solution with a given density profile, while 
we solved the structure of the entire radiative envelope. 

Figure 7 shows theoretical He ZAMS (X = 0) with 
different metal content. The characteristic properties arc 
the same as those of ZAMS in figure 6, but the He ZAMSs 
bend at lower stellar masses. Compared with the helium 
ZAMS obtained by Langer (1989) with the Los Alamos 
opacity, our sequences are close to his curve in the lower 
part (< 5Mq), but bend rightward at less-massive stars 
than does Langer's curve. This is naturally explained by 
the difference in the opacity; the large peak of the OPAL 
opacity causes a core-halo structure in less massive stars, 
which causes bending of the main-sequence. 

This figure also depicts the position of galactic Wolf- 
Rayet stars obtained from a non-LTE radiative trans- 
fer calculation for a helium and nitrogen atmosphere 
(Hamann, Koesterke 1998). Most of the stars distribute 
in the rightside of the curve of Z = 0.02. The open tri- 
angles denote a star in which hydrogen is not detectable 
in the atmosphere. If all of these stars have the same 
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Fig. 7.. Same as those in figure 6, but for helium ZAMS (X = 0). 
The stellar mass at the small dots on the curves are, from lower 
to upper, 5, 8, 20, 30, 40, and 50 Mq for the Z = 0.004 se- 
quence, and 5, 8, 40, and 50 Mq for Z = 0.05. The mass 
of the uppermost point of each curve is 80, 130, 220, and 20 
Mq from left to right. The two horizontal dashed curves con- 
nect the points of 20 Mq (lower) and of 50 Mq on each track. 
The triangles denote the galactic WN stars; the open trian- 
gles represent stars in which hydrogen is not detectable in their 
atmosphere, and the filled triangles with hydrogen (Hamann, 
Koesterke 1998). 

metallicity as that of the sun, Z ~ 0.02, their rightward 
deviation from the curve of Z = 0.02 can be interpreted 
based on evolutional effects. Langer (1989) has followed 
the evolution of helium stars by decreasing the helium 
content of the convective core, and showed that massive 
stars (M > 15M ) leave the main-sequence redward in 
the H-R diagram, while less massive stars (M < 10M Q ) 
move leftward. Considering the difference in the opacity, 
i.e., the large peak in the OPAL opacity is not present 
in the Los Alamos opacity, we can naturally expect that 
stars of mass > 10M Q move rightward in the H-R dia- 
gram. 

The distribution of filled triangles, that denotes stars 
with hydrogen in the atmosphere, seems to be weighted 
to the lower temperature side, which is explained by the 
difference in the opacity. When a helium core-burning 
star has a hydrogen-rich atmosphere, the star intends 
to have a more extended radiative envelope due to an 
increase of the opacity in its surface region. 

A few stars locate to the leftside of the curve of Z = 
0.02, (e.g., a star at logT cff = 4.95 and log L/L Q = 6). 
This may be attributed in part to the difference in the 
definition of the temperature between ours and Hamann 
and Koesterke, in which their temperature, depicted 
here, is defined as that of the bottom of the expand- 
ing envelope, where the optical depth is 20. Therefore, 
the temperature that should be compared with our curve 
may be slightly smaller than these values. In the pres- 
ence of a wind mass-loss, however, a rigid definition of 
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the surface temperature will be very complicated and dif- 
ficult to compare with ours; such detailed examinations 
are beyond the scope of the present paper. Within such 
ambiguities, we therefore conclude that the bending of 
He ZAMS is consistent with the observed data of WN 
stars. 

5. Conclusions 

1. We present the internal structure of chemically ho- 
mogeneous stars with central nuclear burning. Massive 
stars develop a core-halo structure (> 100M Q for MS 
star, and > 15M for helium MS star) with a very ex- 
tended radiative envelope, whereas less massive stars do 
not. This structure change is caused by the large peak 
of the OPAL opacity, which is more prominent in more 
massive metal-rich stars. 

2. The upper part of the main-sequence curves red- 
ward in the H-R diagram, because the core-halo struc- 
ture develops well in more massive stars and the sur- 
face temperature decreases along with an increase in the 
stellar mass. The critical mass of main-sequence turn- 
ing to rightward is 2O(W for Z = 0.004, 100Af Q for 
Z = 0.02, and 70M Q for Z = 0.05 for ZAMS, and 30M Q 
for Z = 0.004, and 15M© for Z = 0.02 for He ZAMS. 
Because the Population I main sequence more strongly 
bends than does Population II, very young massive main- 
sequence stars are expected to distribute in the order of 
the metal content in the upper H-R diagram. 

3. The distribution of observed very massive stars in 
the H-R diagram is consistent with our theoretical main- 
scqucnccs. A ZAMS star of extreme Population I has a 
very low T e g despite its young age, which possibly leads 
to a misclassification as an evolved star. The large radius, 
owing to the core-halo structure of very young MS star, 
must be distinguished from the redward excursion from 
evolutionary effects by observational information on the 
chemical composition. 
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